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At present several groups are analyzing astro-physical data in a search for time- variation of the fine 
structure constant a. Here we discuss how to exclude systematic effects caused by the changes in the 
isotope abundances during the evolution of the universe. We suggest using particular combinations 
of the transition frequencies of O II, Al II, Al III, Si II and Mn II as anchors, which are insensitive 
to a-variation and to changes in isotope abundances. These anchors can be used to determine the 
cosmological redshift. Then, one can use other combinations of frequencies as probes for the time- 
variation of a and another set as probes for the isotopic abundances. In this way it is possible to 
eliminate one source of the systematic errors in the search for a-variation and get new information 
about evolution of the isotopes. On the level of accuracy that has already been reached in the search 
for a-variation it is possible to see ~10% changes in isotope abundances of Si and Ni. 

PACS numbers: 06.20.Jr, 31.30.Gs, 98.80.Ft 



Introduction 

If, during previous stages of the evolution of the uni- 
verse, the fine structure constant a was different from 
its present value ao, we can expect small deviations of 
atomic frequencies from their laboratory values : 



Wlab + qx + 0(x ) 
cj lab + 2qAa/a, 



(la) 
(lb) 



(a/a ) 



1, 



where the parameter q rapidly grows with the nuclear 
charge Z. At least three groups are looking for such 
deviations in the spectra of distant quasars 0,0, EL H, El- 
These groups use different data and their results are not 
fully consistent with each other. In |2j,|3|,|4j a nonzero co- 
variation is reported, while 0, agree with the constant 
a: 



Aa 



-0.57 ±0.11 
1CT 5 x { -0.04 ±0.19 ±0.27. 
-0.06 ±0.06 



sys 
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The frequency shifts which correspond to Eq. J2J are 
of the order of magnitude of the typical isotope shifts. 
Therefore, any changes in isotope abundances can cause 
systematic errors in the search for a- variation (see, for ex- 
ample, discussions in 0,0,11113 and references therein). 



Here we use recent calculations of the isotope shifts 
[ToL ITU Il2l| to reanalyse the role of the isotope evolu- 
tion in the search for a-variation. We suggest how to 
exclude the possible systematic error caused by isotope 
shift from the search for a-variation. We also show that 
the spectroscopic data from Eq. J5J can be used to get 
information about evolution of isotope abundances of Si 
and Ni. 

The shift in frequency of any transition in an isotope 
with mass number A' with respect to an isotope with 
mass number A can be expressed as: 
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where (r 2 ) is the mean square nuclear radius. The mass 
shift constant /cms accounts for the normal mass shift 
and for the specific mass shift [l^: &ms = ^nms + ^sms* 
where 



&NMS 



10 

1823' 



(4) 
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and the value 1823 refers to the ratio of the atomic mass 
unit (amu) to the electron mass. The values of the con- 
stants fcsMS for the specific mass shift and F for the vol- 
ume shift in J21 depend on the details of the atomic struc- 
ture (see ll^lnll f or details). Note that we use the sign 
convention |10L ll T| with &nms < 0, which differs from the 
convention used in |l3j and some other publications. 

Present accuracy of the theory for &sms for atoms with 
more than one valence electron is not always sufficient for 
our purpose. However, it is possible to improve the the- 
ory by applying techniques suggested in 0, 0] . Note 
that calculations for the ions with one electron above 
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closed shells are already sufficiently accurate [lJJ, HB> ■ 
Below we ignore the volume shift of Eq. J3J). For light 
elements it is much smaller than the mass shift, but be- 
comes more important for such elements as Ni, Zn, and 
Ge. Still, its size is comparable to the accuracy to which 
we currently know the isotope shift parameters kyis- 

Isotopic effects 



In this section we discuss the implications of recent cal- 
culations of isotope shifts [TfllTlLll2T | for the astrophysical 
search for a- variation. Table [I] presents calculated values 
of the parameters q and ^Msfor the transitions used in 
astrophysical surveys 0, U S IE El- in order to esti- 
mate possible isotope shifts in astrophysics we calculate 
isotope shifts Aa->i,2 for two simple assumptions: 

• Auj i corresponds to the case when the abundance 
of the leading isotope A is reduced by 10%, while 
the abundance of the next to leading isotope A' is 
increased by the same amount. That means that 
for Si and Cr A' = A + 1; for Mg, Ni, and Zn 
A' = A + 2; for Ti, Fe, and Ge A' = A - 2. Then 
Aid! = 0.1 x Auj a '' A . 

• A0J2 corresponds to the substitution of the terres- 
trial abundances with the single leading isotope. 
For O, Al, and Mn that means zero shifts as these 
elements have no isotopes. 

The shifts AcJi and Atd2 are usually of the opposite 
sign, and we use them to estimate the range of the possi- 
ble isotopic shifts for distant absorber. In order to com- 
pare isotopic effects with the expected effects from co- 
variation, Table^also presents the shifts Aid a which cor- 
respond to Aa/a = -0.57 x 10~ 5 4]. Note that our first 
assumption of a 10% change in abundances is a rather 
conservative one. Sometimes, much more dramatic sce- 
narios of the isotope evolution are considered (see, for 
example Q and the references therein). For these sce- 
narios our conclusions will only be enhanced. 

Table [I] shows that the shifts Acdi^ are indeed of the 
same order of magnitude as the shifts observed in 
Moreover, in some cases there is a strong correlation be- 
tween Auj a and Au>i t 2- That confirms earlier suggestions 
that isotopic effects may lead to significant systematic 
errors in the search for a-variation. 

Let us, for example, consider Fe II. Here isotope shifts 
are relatively small. On the other hand, coefficients q 
and kyis are strongly correlated, and it is impossible to 
disentangle these two effects. For the lines used in |5J, 
the isotope shifts which correspond to our two limiting 
cases would imitate the following a-variation: 



Aa 

a 



= 10 



-5 



-0.14 
-0.07 



case 1, 
case 2. 



(5) 



These values lie within the error bars given by Quast 
et al., but will not allow one to improve their results sig- 
nificantly, unless isotope effects are accounted for. That, 



in turn, will require one to include other lines in the anal- 
ysis. 

Anchors and probes 

When spectra of distant astrophysical objects are ob- 
served, one needs to exclude the cosmological redshift. 
That can be done by looking at the frequency ratios for 
different transitions, i.e. by comparing (u>i/ 'id fc) astro with 
(idi/u)k)\ab- When forming these ratios it is convenient to 
use "anchor" lines, which are not sensitive to the varia- 
tion of the parameter of interest. 

As we saw above, there can be two effects of the same 
order of magnitude, i.e. possible a- variation and variation 
of isotope abundances. Therefore, we need anchors which 
are insensitive to both effects. In Table[|]there is only one 
line of Si II, which may be considered suitable as anchor: 



-4, 
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65495 cm- 



Si ii 



(G) 



It is preferable to have more than one anchor. How- 
ever, all other lines have either relatively large q, or /cmSj 
or both. Aluminum, oxygen and manganese ions have 
only one stable isotope each (leading isotope abundance 
for 16 is greater than 99%). In order to exclude the 
a-dependence of transition frequencies we can use the 
following combinations of 16 O and 27 A1: 
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'aiii 
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■td(3s 1/2 - 
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Aiv 


= 0.68 


• w(3si/ 2 - 
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'aiiii 
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= 19465 cm" 



For 55 Mn we can form several anchors, such as: 



A v = 0.67 • d( 7 S 3 -> 7 P 4 °[3d 5 4p]) MnII (10) 
+ 0.33 • td( 7 S 3 -> 7 P 4 °[3d 4 4s4p]) MnII = 53699 cm" 1 . 

Numerical factors in Q — (|10|l are formed from the q- 
factors for two transitions and are normalized to unity. 
In the leading order An — Ay do not depend on a and 
can be used as anchors. 

Other elements in Table[I]have more than one isotope. 
Strong correlation between the factors q and &ms for Fe II 
allows us to form several combinations of frequencies suit- 
able as anchors. It could be possible to make anchors 
from practically any pair of lines in Fe II which include 
the line 62172 cm -1 , for example: 

Avi = 0.5- W ( 6 C 9/ 2-> 6 ^ /2 ) Fci] 

+ 0.5-w( 6 £> 9/2 -> 6 P? /9 ) = 50316 cm- 1 . 

V ' '"/Foil 
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TABLE I: Results of the calculations of the a- variation parameters q [Til llH . |2(1 l2l| and the isotope shift parameters &ms 
[TflL for the transitions used in astrophysical surveys. The error bars are given in parentheses where available. Column 
Aaj a gives the shifts that correspond to the a- variation observed in p|. Column Ao»i shows the line shift caused by a 10% 
abundance transfer from the leading isotope to the next to leading (by abundance). Similarly, column /S.L02 presents the line 
shift for the case when natural abundance is substituted by a single leading isotope. 



Ion 


Transition 


LOO 


q 


&MS Aio a Auji Alo2 References 






(cm" 1 ) 


(cm" 1 ) 


(cm -1 - amu) (lO^cm -1 ) 



O II 4 S% /2 {2s 2 2p 4 ] -> 4 P 5/2 [2s2p 4 ] 119873 

-»• 4 P 3/2 [2s2p 4 ] 120000 

-> 4 P 1/2 [2s2p 4 ] 120083 

Mgl ^ofis 2 ] -f 1 P?[3s3p] " 35051 

Mg II 3si /2 -» 3pi /2 35669 

-> 3p 3 / 2 35761 

Al II l S [3s 2 ] -» 3 P [3s3p] 37393 

-» 3 Pi[3s3p] 37454 

Al III 3s 1/2 -» 3p 1/2 53682 

-» 3p 3/2 53920 

Sill 2 P° /2 [3s 2 3p] -» 2 L> 3/2 [3s3p 2 ] 55309 

2 5' 1/2 [2s 2 4s] 65495 

Si IV 3s 1/2 -> 3pi /2 71290 

-> 3p 3/2 71750 

Ti II 4 P 3/2 [3d 2 4s] -> 4 F 3 ° /2 [3d 2 4p] 30837 

-> 4 F 5 ° /2 [3d 2 4p] 30959 

CrII e 5 5/2 [3d 5 ] -> 6 P 3 ° /2 [3d 4 4p] 48399 

-> 6 P 5 ° /2 [3d 4 4p] 48491 

-> 6 P 7 ° /2 [3d 4 4p] 48632 

Mn II 7 5 3 [3d 5 4s] -> 7 P 2 °[3d 5 4p] 38366 

-> 7 P 3 [3d 5 4p] 38543 

-> 7 P 4 °[3d 5 4p] 38807 

-> 7 P 2 °[3d 4 4s4p] 83255 

-> 7 P 3 °[3d 4 4s4p] 83376 

-> 7 P 4 [3d 4 4s4p] 83529 

Fell 6 P/ 9/2 [3d 6 4s] -> 6 L»g /2 [3d 6 4p] 6 38459 

-> 6 D? /2 [3d 6 4p] 6 38660 

-> 6 P° 1/2 [3d 6 4p] 6 41968 

-> 6 P 9 ° /2 [3d 6 4p] b 42115 

-> 6 P 7 ° /2 [3d 6 4p] 6 42658 

-> 4 F 7 ° /2 [3d 6 4p] 62066 
-> 6 P 7/2 [3d 5 4s4p] 6 62172 

Nill 2 L> 5/2 [3d 9 ] -> 2 P 7 ° /2 [3d 8 4p] 57080 

-> 2 Z>5 /2 [3d 8 4p] 57420 

-> 2 P 5 ° /2 [3d 8 4p] 58493 

Zn II 4s 1/2 -> 4p 1/2 48481 

-> 4p 3/2 49355 

Ge II 4p 1/2 -> 5si /2 62403 
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[10, .19] 


2490 


(25) 


-69.3 (2.3) 
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"For this transition we give experimental value for k^s from |22| . 
'The lines of Fe II used by Quast et al. 0. 
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Now, when we have several anchors to exclude the cos- 
mological redshift, we can look for the combinations of 
frequencies which are sensitive to only the mass shift or 
to the variation of a. These combinations can be used 
as "probes" of isotope abundances and a-variation. For 
example, the combination 



TABLE II: a-variation probes, insensitive to isotope shifts. 



Pi = 0.64-w(3a 1/2 -» 3p 1/2 ) MgH 

- 0.36 • w(3« 1/2 -» 3p 3/2 ) Mgn = 9740 cm" 



(11) 



is insensitive to a-variation and can serve as a probe of 
the abundances of Mg isotopes in distant absorbers. 

Note that even though g-factors for Mg are relatively 
small, it is still necessary to use Pj as a probe, rather 
than individual lines. Since a-variation on the scale of 
Aa/a ~ —0.6 x 10 -5 is currently not excluded, we can ex- 
pect the shift <5w(3s 1 / 2 — + 3p 3 / 2 ) ~ —0.002 cm -1 , which 
is comparable to the current sensitivity of observations 
to the frequency shifts. 

By analogy with Mg II it is possible to form combina- 
tions for the fine structure doublets in Si IV and Zn II: 



ii 



in 



= 0.71 • w(3s 1/2 -> 3pi /2 ) 

- 0.29 • w(3si/ 2 -> 3p 3/2 } 
= 0.61 • w(4s 1/2 -> 4p 1/2 ) 

- 0.39-w(4s 1/2 ->4p 3/2 ) 



Si IV 



Si IV 



= 30381 cm' 



Znll 



(12) 



(13) 



Znll 



= 10423 cnT 



For most other transitions in Table H] there is strong 
correlation between parameters q and fcMSj an d we can- 
not exclude the dependence on a. The only exception is 
Ni II, where the transition with frequency 58493 cm -1 is 
insensitive to a-variation and two other transitions have 
different parameters q. That gives us two more probes of 
isotope abundances: 



Piv 
Pv 



u( 2 D 5/2 -► 2 F 5 ° /2 ) = 58493 cm" 1 , (14) 

V / Ni II 

0.67-o;( a D 6/2 -> a F?' /3 ) (15) 



Nill 



0.33-w( 2 L> 5/2 -> 2 D° /2 ) = 18952 cm- 1 . 



Nill 



The same method can be applied for building combi- 
nations that are insensitive to isotope shift and can serve 
as probes for the variation of a. Above we suggested that 
several combinations of spectral lines of O, Al, and Mn 
ions can be used as anchors. These lines themselves can 
also be introduced as a-dependent probes. Indeed, these 
elements have only one stable isotope each and there- 
fore have no isotope shifts. Here, manganese is the best 
probe, as it has the biggest absolute value of q. Ge II is 
also a good probe because of a very low value of /cms of 
the only transition 4p 1 / 2 ~* 5si/2> listed in Tabled 

Other ions from Table |U are sensitive to isotope shifts. 
Still it is possible to make several combinations to ex- 
clude this dependence. Mg II, Si II, and Ti II have rela- 
tively low values of q and would not make good probes. 
For remaining ions fcMS is practically the same for all 
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Al III 
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w(si/ 2 -> 3p 3/2 ) 
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-1/2 ■ uj(3s 1/2 -> 3pi/ 2 ) 

1 In In ro lb / .IN 
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38807 
83255 
83376 
83529 
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-3033 
-2825 
-2556 
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0)( 7 S , 3 

uj( 7 S 3 



7 P 3 [3d 5 4p] 
■ 7 P 4 [3d 5 4p] 
7 P 2 [3d 4 4s4p] 
7 P 3 [3d 4 4s4p] 
7 P 4 [3ri 4 4s4p] 
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-350 


1/2. 


wC-Db/2 


^ ^5/2) 
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-1/2 
1/2 
-1/2. 


w( 2 I>5/2 

oj( 2 D 5/ 2 

w( 2 A>/2 


t 7/2) 
-* ^5/2) 
u B/2) 


Zn II 
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453 


1/2 
-1/2 


■ w(4si/ 2 

■ w(4si/ 2 


-» 4p 3/2 ) 
-» 4pi /2 ) 


Ge II 


62403 


607 




w(4pi /2 


-> 5s 1/2 ) 



transitions of interest. Hence, we can form a probe by 
taking the difference of any two frequencies. The results 
are given in Tabled where we introduced normalization 
factors as above. 

For Fe II and Cr II it is impossible to exclude isotope 
dependence because of the correlations between q and 
/cms- In both cases isotopic effects are relatively small. 
That allows one to use these ions to test a-variation on 
the level reported in [j] , but not on the level of the result 



Discussion of sensitivity 

Let us use Eq. J2J and Table |U to estimate typical sen- 
sitivity to the frequency shifts in modern astrophysical 
studies. The paper |f| is based only on the analysis of the 
lines of Fe II. From Table[Q we see that for all transitions 
in this ion \q\ « 1300 cm -1 . Eq. Ijlbf) then gives statisti- 
cal accuracy Alo w 0.005 cm -1 . Two other groups 0,0 
use more ions, but Fe II is present in both samples and 
we can again use the same value of |g| as typical. That 
leads to the estimates of the accuracy which are roughly 
2 Q and 3 [(| times better than above. This is statistical 
accuracy, which is achieved by averaging results for many 
transitions, so we will use the following estimate of the 
achievable accuracy for the frequency shifts in astrophys- 
ical observations: 

Aoj w 0.005 cm" 1 . (16) 
Now we estimate frequency shifts for the probes Pjr - 
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Py and compare them with current sensitivity flS} . To 
do this we need to specify the way in which isotope abun- 
dances are changed. We assume that the intensity of the 
line of the leading isotope is transferred to the next to 
leading one. If there are two comparable weaker lines, we 
assume that they are equally increased. Within this sim- 
ple model we calculate the shifts of the centers of gravity 
of the lines and the resulting shift of the probe. Com- 
parison of this shift to the modern frequency sensitivity 
(|To)) gives us the sensitivity to the isotope abundances. 

Mg has three stable isotopes A = 24, 25, and 26 with 
modern abundances 79 : 10 : 11. Suppose that in the 
early universe it was (79 — x) : (10 + x/2) : (11 + x/2). 
Using factors /cms from Table Hand ignoring the volume 
shift and saturation effects, we can estimate correspond- 
ing shift of the center of gravity of the probe l|ll|) : 



6P T = -0.27 



3fcMS^ 

200A 2 



0.00023 x cm" 



(17) 



For Si and Zn there are also two comparable admix- 
tures to the leading isotope. If we assume (92 — a;) : 
(5 + x/2) : (3 + x/2) abundances for Si and (49 - y) : : 
(28 + y/2) : 4 : (19 + y/2) : : 1 abundances for Zn, we 
get: 



SPn 
SPm 



(,.43^M^« 0.00074a; cm" 1 . (18) 



-0.22 



200 28 2 
3k MS y 

100 64 2 



0.00011 y cm- 1 . (19) 



Ni has only one dominant admixture to the leading iso- 
tope, so we assume (68 — x) : : (26 + x) : 1 : 4 : : 1 
abundances for Ni. That gives: 



5P 



IV 



2k M s x 
100 58 2 



0.00046 a; cm" 



(20) 



while the probe Py is three times less sensitive. 

Comparing (|17fl - l|20|) with (|16fl we see that current 
sensitivity allows detection of a 22% change in the abun- 
dance of the isotope 24 Mg, a 7% change for the isotope 



28 Si, 45% for 64 Zn, and an 11% change in the abundance 
of 58 Ni. 



Conclusions 

In this paper we analyze the sensitivity of modern as- 
trophysical surveys that search for a- variation to isotopic 
effects and come to the following conclusions: 

• Isotopic frequency shifts are generally likely to be 
of the same order of magnitude as the present ex- 
perimental sensitivity. Therefore, isotope shifts are 
an important source of systematic errors. 

• It is possible to eliminate isotopic effects by tak- 
ing proper combinations of transition frequencies. 
These combinations depend on the calculated mass 
shift coefficients. At present not all of them are 
known with a sufficient accuracy. Corresponding 
theory can be improved by applying the method 
suggested in [l4l Il5| . New experimental data on 
the isotope shifts is also highly desirable. 

• As a by-product of the search for a-variation it is 
possible to get information about isotopic evolution 
of several elements. 



• The accuracy of the surveys, which use only lines 
of Fe II is limited by the isotopic effects on the 
level 8a/ a ~ 10~ 6 . This estimate is based on a 
rather conservative assumptions about possible iso- 
tope evolution. 
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